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Abstract 

U , 

Using the multidimensional dynamic-programming method (MDPM) 
in the four-dimensional deformation space {f3\} with A=2, 4, 35 and 6 
we were able to study evolution of the action integral of the fissioning 
CN ! nucleus 258 Fm. We found the second minimum on the cross-section of 

the action integral for /?2 ~ 1 , what we interpret as the dynamical 
evidence of the bimodal fission in this heavy Fm isotope. 
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PACS numbers: 25.85. Ca, 24.75.+i, 02.60. +y 
^ : 1 Introduction 



Experimentally known values of the spontaneous fission half-lives T s f of nine 
even-even Fm isotopes (N = 142, 144, 158) form approximately two sides 
of an acute-angled triangle with a vertex in N = 152. The rapid changes in T s f 
on both sides of 252 Fm are particularly dramatic for the heavier Fm isotopes, 
where T s f descends by about ten orders of magnitude when one passes from 
254 Fm to 258 Fm. 

^ ■ Nuclei in the vicinity of 258 Fm reveal another two peculiarities in spon- 

taneous fission properties. The fission process becomes symmetric with the 
narrow mass distributions (see e.g. review [l|]) and the total-kinetic-energy 
(TKE) distributions are not Gaussian, but instead are best described as a 
sum of two Gaussians. For 258 Fm the TKE distribution has peaks at about 
200 and 235 MeV. It was postulated [0 that for 258 Fm the observed high- 
kinetic-energy peak corresponds to fission through a scission configuration of 
two touching spherical fragments, with the maximum of Coulomb repulsion. 
The low-kinetic-energy peak was thought to correspond to fission through 
a conventional scission configuration of two elongated fragments. 

The investigations focused on a multi-valley structure of potential-energy 
surfaces for nuclei near Fm have been undertaken by a lot of authors. A par- 
ticularly detailed studies were performed by a Los Alamos group |3|, £|, |j| and 
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a Warsaw group || Comparison of used models and results obtained by 
the both groups one may find in || and in a review ||. 

The descriptions of the bimodal fission in 258 Fm proposed by the Los 
Alamos and Warsaw groups are very similar. Roughly speaking in both cases 
authors present two fission path on the potential-energy surfaces. One trajec- 
tory goes into the valley which corresponds to more compact shapes (CS) of 
the nucleus and the other to the valley which corresponds to more elongated 
shapes (ES). However, the both presented descriptions of the bimodal fission 
have one common deficiency — they are totally static, focused only on the 
multi-valley structure of potential-energy surfaces. 

The main goal of this paper is to present a dynamic analysis of the bimodal 
fission in 258 Fm. Thus, we focused our attention on an action-integral of the 
fissioning nucleus in a multi-dimensional deformation space. 

The used model is described in sect. 2, the results and discussion are given 
in sect. 3. 



2 Model 

To calculate the action-integral along a path L(s) in the multi-dimensional 
deformation space {X\} we used the one-dimensional WKB semiclassical ap- 
proximation 

S(L) = J" l [ ^B cS (s)[V(s)-E]j 1 \s, (1) 

where an effective inertia associated with the fission motion along the path 
L(s) is 

A, /i 

In above equations ds defines the element of the path length in the {X\} 
space. The integration limits S\ and S2 correspond to the entrance and exit 
points of the barrier V(s), determined by a condition V(s) = E, where E is 
the penetration energy of the fissioning nucleus. 

The potential energy V is calculated by the macroscopic-microscopic model 
exactly the same like that used by the Warsaw group 0. For the macro- 
scopic part we used the Yukawa-plus-exponential finite-range model || and 
for microscopic part the Strutinsky shell correction, based on the Woods-Saxon 



single-particle potential with "universal" variant of the parameters [1D|. The 
single-particle potential is extended to involve residual pairing interaction, 
which is treated in the BCS approximation. The inertia tensor Bx x x x > which 
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describes the inertia of the nucleus with respect to change of its shape, is cal- 



culated in the cranking approximation (cf. e.g. fll!]). The penetration energy 
E = V(X®) + 0.5 MeV is defined as a sum of a ground-state energy V(Xy) 
at the equilibrium deformation X® and a zero-point energy in the fission di- 
rection, equal 0.5 MeV. 

Dynamic calculations of the spontaneous fission half-lives T s f are under- 
stood as a quest for a least-action trajectory L min which fulfills a principle 
of least-action ^[^(L)] = 0. To minimize the action integral (1) we used 
the multidimensional dynamic-programming method (MDPM). Originally this 
method was used only for two-dimensional deformation space [pjjj . We ex- 
tended the model up to four degrees of freedom. 

The schematic Fig. 1 demonstrates how our model, the MDPM, works. 
Since the macroscopic-microscopic method is not analytical, it is necessary 
to calculate the potential energy and all components of the inertia tensor 
on a grid in the multidimensional space spreaded by a set of deformation 
parameters {^a}- We select one coordinate X from this set. This coordinate 
(e.g. elongation parameter) is related in a linear way to the fission process. 
In Fig. 1 for each point on a X Q axis the rest of the coordinates {X\_i} are 
represented for simplicity by a "wall-plane" of only two coordinates X\ and 
X 2 . 

To find the least-action trajectory L min between the turning point s\ and 
$2 we proceed as follows. First, from the entrance point to the barrier s\ we 
calculate the action-integrals to all grid points in the nearest "wall" at X = 1. 
In the next step we come to the "wall" at X = 2 and from each grid point 
in this "wall" calculate the action-integrals to all grid points in the "wall" at 
X = 1 (see Fig. 1). The trajectories started from each grid point at X = 2 
passing through all grid points in "wall" at Xq = 1 and terminated in the point 
Si form a bunch of paths. From each such a bunch we choose the path with 
minimal action-integral and keep them in the memory. At the end of this step 
we have the least-action integrals along trajectories which connect the starting 
point Si with all grid points in the "wall" at X = 2. After that we repeat this 
procedure for all grid points at X = 3 and again we obtain all least-action- 
integrals along trajectories starting from point Si with ends at each grid point 
in the "wall" X = 3. We repeat it until we exceed the n-th "wall" ; a last one 
before the exit point from the barrier S2- After that we proceed to the last 
step of our method. We calculate action-integrals between all grid points in 
the last "wall" at Xq = n and the exit point S2] the minimal one among them 
corresponds to the searched trajectory of the least-action-integral L min . 

Our dynamical calculations are performed in the four-dimensional defor- 
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mation space spread by /3-shape parameters appearing in an expansion of a 
nuclear radius in spherical harmonics. So, we have {X\} = {P\}, with A = 
2,4,35 and 6. Our deformation space is very similar to used in ref. 0. How- 
ever, in opposite to ref. @, we collect the reflection-asymmetry parameters 
Ps and /?5 in one parameter j3 35 = (/3 3 , P5—O.8P3), according to the average 
trajectory in a (P3, (3§) plane for nuclei in Fm region. 

The potential energy V and ten components of the symmetric inertia tensor 
Bp x [3 x are calculated microscopically for 

fa = 0. 15(0.05)1.30, 1 

P A = -0.08(0.04)0.36 , (3) 
p 35 = 0.00(0.05)0.25 , 
Pe = -0.12(0.04)0.12 

i.e. in the 24 x12x6x7 = 12096 grid points. In our calculation a quadrupole 
deformation p 2 plays a role of the coordinate X from Fig. 1. 

3 Results and discussion 

Fig. 2 displays the cross-sections of the action-integral (in units of h) of the 
nucleus 258 Fm, obtained according to our MDPM model for the different values 
of parameter p%. To show a three-dimensional structure of these cross-sections 
we reduced the results to the two-dimensional (/?4, ^35) contour plots. This we 
accomplish by plotting at each point a minimal value of action-integral with 
respect to a third coordinate Pq. Thus, in Fig. 2 one can see an evolution of 
the action-integral in the fission process of 258 Fm. 

The two-dimensional contour diagrams show fairly smooth surfaces with 
one minimum. On all contour diagrams a position of the minimum, in a lower- 
left corner, is almost the same. If we join the minima from all contour maps, 
we obtain a dynamical fission trajectory in our four-dimensional deformation 
space. The dynamical fission trajectory has a tendency to be close to a straight 
line. It means a reduction of the effective inertia -B e ff according to eq. (2) and 
leads to the smaller action-integral. This dynamical fission trajectory goes to 
the compact shapes (CS) (cf. e.g. Fig. 6 in ref. J7J). 

In the last contour map in Fig. 2 with the cross-section of the action- 
integral for Pi = 1.05 one can see in an upper-right corner a second local 
minimum corresponding to more elongated shapes (ES). This second local 

1 Here we use the notation x — Xi(Ax)xf which means that x goes from a;, through Xf 
with step Ax. 
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minimum is more distinctly seen in Fig. 3 where we plot a contour diagram of 
the action-integral obtained in the last step of the MDPM, when we calculate 
the action-integral from all grid points on last "wall" to the exit point from 
the barrier s 2 (see Fig. 1). In other words, in Fig. 3 we show a "view" of the 
last "wall" at fa — 1.05 "seen" from the exit point s 2 with coordinates: (fa = 
1.10, fa = 0.05, fa 5 = 0.01, fa = -0.06). The appearance of two minima on the 
cross-section of the action integral, one global corresponding to more compact 
shapes (CS) and another local corresponding to more elongated shapes (ES), 
indicates that in the vicinity of fa = 1.05, before the exit from the barrier, one 
can observe the begining of the two fission modes in 258 Fm . 

Fig. 4 shows the potential energy V calculated for 258 Fm as a function of 
the deformation parameters fa and fa. At each point (fa, fa) we select the 
value of V(fa, fa, (3 35 , fa) minimal with respect to /5 35 and fa, where /3 35 and 
fa are taken only from the grid points, according to (3). On the right part of 
the figure one can see two valleys. Entrances to the valleys of compact and 
elongated scission shapes are denoted by CS and by ES, respectively. 

Fig. 4 also includes paths to fission. A nearly stright, solid line, which 
starts in the first minimum (at fa ~ 0.25, fa ~ 0.0) and goes to the CS valley, 
represents the dynamical fission trajectory leading to compact shapes. From 
this trajectory, close to the exit point from the barrier (at fa ~ 0.95, fa ~ 0.05), 
the fission trajectory leading to the ES valley branches off (a dotted-dashed 
line). 

Only for comparison we display in Fig. 4 a static fission trajectory (i.e. the 
path of minimal potential energy for each fa value), shown as a dashed line. 
This trajectory starts in the first minimum, passes through the first saddle 
point (at fa ~ 0.4, fa ~ 0.1), then through the second minimum (at fa ~ 
0.65, fa « 0.05), through the second saddle point (at fa « 0.85, fa « 0.05), 
and leads to the valley of elongated scission shapes. 

All fission barriers along the static and both dynamical trajectories in Fig. 
4 have two humps, in distinction to ref.J7[ where the static barrier has only 
one hump. This difference, despite the similar macroscopic-microscopic models 
used, seems to come from minimization procedure of the total energy of the 
nucleus in the multidimensional deformation space applayed by the Warsaw 
grup in §]. 

Our results of bimodal fission in 258 Fm are in agreement with those obtained 
by the Los Alamos group ||, |, |j, §], where the splitting of the trajectories 
appears also a little befor the exit point from the barrier. This is in opposite to 
the Warsaw group || |7j description where the fission trajectory splits behind 
the exit point from the barrier. 
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Figures Captions 



Figure 1: Schematic presentation of the MDPM method. In the multidimen- 
sional deformation space {X\} we select the coordinate {Xq} which is related 
in a linear way to the fission process. The points Si and s 2 correspond to 
entrance to the barrier and exit from the barrier, respectively. See text for 
details. 



Figure 2: Evolution of the action-integral in the fission process of 258 Fm. 
Each contour map represents cross-sections of the action-integral, in h units, 
obtained for the different values of parameter /3 2 - 



Figure 3: View of the last contour map from Fig. 2, with cross-section of 
the action-integral at $2 = 1.05, "seen" from the exit from the barrier point 
S2- The letters CS indicate global minimum corresponding to more compact 
shapes, ES indicate second local minimum corresponging to more elongated 
shapes. 



Figure 4: Contour map of the potential energy of 258 Fm, showing paths to fis- 
sion. Solid line represents the dynamical fission trajectory leading to compact 
shapes (CS), the dotted-dashed line denotes branches off second dynamical 
trajectory leading to elongated shapes (ES). The dashed line represents the 
static fission trajectory. 
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